Magnetic-field asymmetry of nonlinear mesoscopic transport 
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We investigate departures of the Onsager relations in the nonlinear regime of electronic transport 
through mesoscopic systems. We show that the nonlinear current- volt age characteristic is not an 
even function of the magnetic field due only to the magnetic-field dependence of the screening 
potential within the conductor. We illustrate this result for two types of conductors: A quantum 
Hall bar with an antidot and a chaotic cavity connected to quantum point contacts. For the chaotic 
cavity we obtain through random matrix theory an asymmetry in the fluctuations of the nonlinear 
conductance that vanishes rapidly with the size of the contacts. 
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Introduction. — The Onsager-Casimir relations Q, Q 
are symmetry conditions for correlation functions. These 
can be cast as friction coefficients by means of the 
fluctuation-dissipation theorem. In electronic transport 
measurements, the information about dissipation is car- 
ried by the conductance of the sample. Microscopic re- 
versibility thus requires that in the presence of a mag- 
netic field B the conductance obeys G a p(B) = Gp a (—B) 
between contacts a and /3. In particular, for a two- 
probe conductor the conductance is an even function of 
magnetic field G(B) = G(—B). Such relations gener- 
ally hold for macroscopic systems near thermodynamic 
equilibrium. What happens when these conditions are 
not fulfilled? First, when transport is phase coherent as 
it occurs in mesoscopic conductors, the conductance is 
not just material specific but also depends on the probe 
configuration. Then, a generalized reciprocity theorem 
may be proved provided current and voltage probes are 
treated on equal footing 0. Second, away from equilib- 
rium (e.g., in nonlinear dc-transport) there are no fun- 
damental reasons why the Onsager relations should hold. 
In particular, in a two terminal conductor, we can expect 
that the nonlinear I—V characteristic is not an even func- 
tion of magnetic field, I(B,V) ^ I(—B,V). Below, we 
demonstrate that such asymmetries can indeed be found 
but importantly only to the extent that the potential 
landscape in the interior of the conductor becomes an 
uneven function of magnetic field Q. 

Nonlinear transport in mesoscopic conductors is a sub- 
ject of growing interest. Rectification effects in asymmet- 
ric micro junctions quantum point contacts p, and 
switching and gain in three-terminal ballistic branches 
have been observed. For samples with special spatial 
symmetries, a set of symmetry relations which hold in the 
nonlinear regime have been discussed and experimentally 
verified by Lofgren et al. || Conductance fluctuations 
away from equilibrium have long been of theoretical 
and experimental [ToL Hlj concern. Recently, Zumbiihl, 
Marcus et al. 0] started experiments to investigate the 
magnetic-field symmetry of fluctuations away from equi- 
librium. 



A correct treatment of this problem requires a self- 
consistent discussion of the screening potential within 
the conductor. A scattering picture of weakly nonlinear 
transport was provided by Christen and Biittiker 0, ^| . 
The work of Sheng et al. [14j exemplifies a computational 
implementation. We briefly now discuss the essential el- 
ements of this approach. 

Consider a mesoscopic conductor penetrated by a mag- 
netic field B perpendicular to the plane of the sample 
and connected to a — 1 , . . . , M reservoirs and gates. 
Transport is described by the scattering matrix s a p . For 
purely elastic transport, the scattering matrix is a func- 
tion of the energy E of the carriers and it is a func- 
tional of the potential landscape U(r) of the conduc- 
tor, s a p[E;U(r)]. For linear transport, the scattering 
matrix is evaluated at the equilibrium potential U cq (r) 
which is an even function of B. As a consequence, 
s a fj{B) = Sfj a {—B) at equilibrium. Away from equi- 
librium the potential depends on the voltages V a (mea- 
sured from an equilibrium chemical potential ^o) ap- 
plied to the leads and the nearby gates. We can write 
U(f) = U cq (r) + J2 a u a V a + 0{V 2 ). Here, the character- 
istic potentials (CP's) u 7 (r) — [dU (f) / dV 7 ] eq relate the 
variation of U(r) in the sample to a voltage shift in the 
contact 7 0. We expand the current through lead a in 
powers of the voltage shifts up to second order: 
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The linear conductances are expressed by the well 
known formula G a p = (e 2 /h) J dE[—dEf(E — /j,q — 
Va)]A a /3(E; {V a } = 0), where f(E) is the Fermi function 
and A a p(E;{V a }) = Tr[l a p5 a p - s' a/3 s a p]. The second- 
order nonlinearity G a p-y contains information about the 
charge response of the system. Reference [13 ] finds: 
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2u 7 (r)]. (2) 



Here the nonequilibrium state is described by the CP's, 
up, which arise as a consequence of screening of the addi- 
tional bare charge injected from contact f3. This charge 
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density of states (DOS) is the injectivity Q of contact /?, 

dfl ^) ^ / , 8s al3 \ 

dE ~ 2m ^ { S ^e5U(f)J- ( ) 

On the other hand, the additional current in contact (3 
due to a variation of the screening potential at point r, 
is given by the emissivity into (3, dnp{r)/dE. Due to mi- 
croreversibility, dnp(B,r)/dE = dn^(~B,r)/dE. How- 
ever, neither the injectivity nor the emissivity alone are 
invariant under i?-reversal. As a result, the CP's are not 
even functions of B. Thus, quite generally, even in a two- 
terminal setup the second-order contribution G a p~ ( to the 
I—V characteristic is not even in the field. We define the 
magnetic- field asymmetry for such a setup: 

$ = ^[Gm(B)-G m (-5)]. (4) 

We emphasize that the asymmetry $ ^ is generated in 
Eq. @ only through the asymmetry in the electric po- 
tential: if the potential is even in B, Eq. (J2J) predicts for 
a two-terminal conductor a current that is even. Thus, a 
self-consistent description of charge redistribution is cru- 
cial. 

Our purpose is to elucidate this general result with the 
help of a simple but instructive example (a quantum Hall 
bar with an antidot) and to provide a prediction of the 
size of $ for a generic conductor, a chaotic cavity. 

Quantum Hall bar. — We consider a conductor in the 
quantum Hall regime as depicted in the left panel of 
Fig. ^ Edge states are symbolically indicated by arrows 
along the upper and lower edge of the sample. For sim- 
plicity, we assume that B is so strong that only the lowest 
Landau level is occupied (filling factor v — 1). Backscat- 
tering is achieved by producing with gates a potential 
hill, thus forming an antidot |l5| . The antidot behaves 
effectively like a quantum impurity with a Breit-Wigner 
resonance at Eq + eUd coupled to the edge states via 




FIG. 1: Left panel: quantum Hall conductor attached to two 
reservoirs with an antidot connecting two edge states. Right 
panel: Schematic spatial variation of the screening potential 
inside the bar for opposite polarities of the magnetic field. 



hybridization widths Ti and L2 due to tunneling. Both 
broadenings are taken to be energy independent. Ud is 
the potential at the antidot. 

Let us first provide an intuitive picture of the nonlinear 
transport. Application of a voltage V\ > V2, leads to an 
excess charge on the upper edge, and through screen- 
ing to a corresponding deficit of charge on the lower 
edge. Hence, a Hall potential r) is established. 

The charges on the edge state lifts the resonant energy 
to a value Eq + eUd(B) (upper right panel of Fig. 
If we now reverse the field, it is the lower edge state 
that is charged, and the upper edge state that is low- 
ered in energy through screening, generating an electric 
field opposite to that of Uh(B). Then, the resonant en- 
ergy will now be at E + eUd{—B) (lower right panel of 
Fig.^J. Clearly, elld(B) and eUd(-B) will in general not 
be equal, except if stringent symmetry conditions are ful- 
filled. The symmetry is broken either through scattering 
asymmetry, i.e., if transmission into the antidot is not 
symmetric Li ^ T2, or through electrical asymmetry if 
the charges on the upper edge couple more strongly to 
the antidot than the charges on the lower edge. 

We now support this picture with analytical calcula- 
tions. Evaluation of the exact local potential distribution 
can typically be achieved only computationally. To sim- 
plify the problem, we divide the conductor into five re- 
gions tti with i — 1, . . . , 5 as indicated in Fig.^ S7s = fid 
is the region of the edge state circling the antidot. In each 
region the potential is taken to be constant. Interaction 
between charges in different regions is described by a geo- 
metric capacitance matrix dj . Such a discrete potential 
model captures the essential physics 

HG3. Thus, the CP 
in region i is Ui 7 — dUi/dVj and the injectivity of lead a 
into region i is Di a = J„ d 3 r dfi a (r)/dE. When the en- 
ergy is close to E + eUd an electron is reflected through 
the antidot (e.g., from fli to ^3) with a probability R = 
1-T = rir 2 /|A| 2 , where A = ^ - E - eU 5 + iT/2 with 
r = Ti +T-2 as the total linewidth. Using the correspond- 
ing scattering matrix we find the injectivities into the dif- 
ferent regions. For the case B > (Fig.Q] left panel) the 
DOS associated to carriers at injected by lead 1 reads 
D31 = D3R, where Di is the total (independent) DOS 
of the upper edge state in region i. Similarly, the remain- 
ing injectivities of contact 1 are D\\ = D-y, D21 = D2T, 
D41 = and D 5 i — e 2 Ti/2n\ A| 2 . One proceeds likewise 
to find the injectivities of contact 2, D^- The charge 
in region i can be expressed in two ways: 

9i = e 2 J2 Ac* {V a -Ui) = J2 Qj Uj . (5) 

a j 

First, qi is the bare injected charge due to the voltage 
applied to the contact and the screening charge induced 
by the internal potential Ui. Second, qi is the charge 
permitted by the Coulomb interaction, where Cy is the 
geometrical capacitance matrix whose indices run over 
all (five) regions considered in Fig. 2] Equation (J3J) al- 
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lows us to determine the potentials J7, as a function of 
the applied voltages. We take equal DOS for all regions 
Di = D. We consider separately the case of (A) an elec- 
trically symmetric sample that is asymmetric only in the 
scattering properties and (B) the case of a sample with 
symmetric scattering that is asymmetric only electrically. 

For case (A) we assume a capacitance matrix with 
equal capacitances C between the edge states and the 
antidot. Transmission from the antidot to the upper and 
lower edge state is asymmetric 77 = (ri — T2)/T. Us- 
ing Eq. J5J and Eq. (0 we evaluate the second-order 
nonlinear conductance and insert it in Eq. (@}. Taking 
into account that the transmission (reflection) probabili- 
ties depend only on the potential at the antidot, and that 
ST/eSUd — —dT / <1E\e f at zero temperature, we find 



$ = - 



h dE 



r]R(C + e 2 D) 



2nCDT + R{C + e 2 D) 



0(V 3 )- (6) 



We observe that in the charge neutral limit (C = 0) 
the magnetic-field asymmetry $ is independent of R (to 
leading order in 77). In the nonintcracting limit (C — > 00), 
the asymmetry is proportional to R/(R + 2ttDT). 

Next consider case (B), a sample with symmetric scat- 
tering properties (77 = 0) but being electrically asymmet- 
ric. Such a case arises if the tunnel barrier separating the 
antidot from the top edge state is not very high in energy 
but wide whereas the barrier separating the antidot from 
the lower edge state is high in energy but narrow such 
that transmission through both is equal. However, the 
capacitance to the upper (1+£)C and lower (1 — £)C edge 
states will differ, £ being a dimcnsionlcss parameter. As 
above, we calculate the -B-dependent second-order con- 
ductance to lowest order in £ to find: 



~ ~~h dE 



TT^e 2 D 2 CTT 



cq 



(C + e 2 D)[2nCDT + R{C + e 2 D)} 



(7) 

Interestingly, the magnetic-field asymmetry is propor- 
tional to the transmission probability T, unlike case (A). 
Thus, the transmission of the impurity would act exper- 
imentally as the indicator of the physical mechanism be- 
hind the resulting field asymmetry. 

Our discussion demonstrates that either asymmetric 
scattering or an electrical asymmetry generates already 
to second order in voltage a deviation from the On- 
sager relations that hold in the linear regime. For strong 
backscattering, this deviation is due mainly to a scat- 
tering induced asymmetry whereas for weakly coupled 
impurities the electric asymmetry dominates. 

Chaotic cavity. — It is important to find out whether 
such symmetry breaking is only relevant for few-channel 
problems or whether in fact symmetry breaking is also 
observable in the many-channel limit. This is the motiva- 
tion to investigate now a chaotic cavity, which is a metal- 
lic quantum dot whose classical analog displays chaotic 



dynamics. For clean samples the transport within the 
cavity is ballistic and its corresponding statistics is well 
described by random matrix theory [lfj]. Open cavities 
have been extensively studied both theoretically 0] and 
experimentally Most of these works are restricted 

to the linear conductance regime (see, however, Ref. ||). 
The chaotic cavity is coupled to reservoirs (a = 1,2) 
through quantum point contacts with Ni and N 2 prop- 
agating channels (left inset of Fig.JSJ. On the ensemble 
average such a cavity exhibits simply a linear I— V char- 
acteristic with a conductance G = e 2 NiN2/hN where 
N = Ni + N 2 . Nonlinearities arise due to quantum fluc- 
tuations with an energy scale equal to the Thouless en- 
ergy Et = N8 with 6 the mean level spacing. 

Since a cavity is effectively zero-dimensional due to its 
isotropic properties, we take into account screening with 
a single potential U. Surrounding gates are coupled ca- 
pacitively with a gate voltage V g and a geometric capaci- 
tance C. In response to a shift of the contact voltages V a 
a noncquilibrium charge builds up in the cavity which de- 
pends on the injecting contact. The two equations which 
determine the excess charge on the cavity are 



= e 2 D a (V a -U) = C(U - V g ) . 



(8) 



The 



where, as before, D a is the injectivity of lead a. 
total DOS of the cavity is D = ^2 a D a . 

As we have now a three-terminal problem with a gate 
voltage V3 = V g , we consider the case where V3 = V 2 - 
FromEq. Q we get h = G 11 {V 1 -V 2 ) + G 111 (V 1 -V 2 ) 2 + 
0(V 3 ). Using the WKB approximation, we can again 
replace derivatives with regard to potentials by energy 
derivatives |l3( and find 



e 3 dT 

Glll -~TdE 



(l-2ui). 



(9) 
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FIG. 2: (Color online). Normalized probability densities for 
the magnetic-field asymmetry $ of a chaotic cavity connected 
to two reservoirs (left inset). Right inset: fluctuations of $ as 
a function of the total number of modes N. The dashed line 
is the analytical prediction for N 3> 1. 
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For the chaotic cavity the ensemble average (Gin) = 
vanishes, an asymmetry can develop only due to quantum 
fluctuations. We now show that var($) is nonzero. 

For Ni 3> 1 and N 2 3> 1 within random matrix the- 
ory the product in the right hand side of Eq. JHJ can 
be decoupled 0. Thus, we disregard correlations be- 
tween dT/dE and u\. The unscreened nonlinear con- 
ductance — (e 3 /h)dT/dE\ eq changes sign randomly on 
the ensemble so that its average is zero 0|. In the 
unitary ensemble (magnetic flux through the cavity of 
the order of 1 quantum) we find for the fluctuations 
vax(dT/dE) = 8ir 2 NfN$/N 6 5 2 . Furthermore, we ne- 
glect the small fluctuations hi the electrochemical 
capacitance 1/G M = 1/C + S/e 2 . Then, from Eq. (JHJ and 
Eq. JHJl we find that the fluctuations of u\ are determined 
by the fluctuations of the asymmetric part of the injec- 
tivity. Correlations of injectivities have been investigated 
in Ref . [2(J ■ Using these results in Eq. we find 

The fluctuations go as 1/N 4 for N ^> 1, thus vanish- 
ing quickly with increasing number of channels and are 
maximal for perfect screening G M = C. 

For small Ni and N 2 a full analytical calculation of 
var($) is involved since one has to evaluate high-order 
correlations between D a , D, and dT/dE. Yet, we can 
substantiate Eq. pO(l with numerical calculations. The 
results for C = and N% = N 2 = N/2 are shown in 
Fig- El I n these simulations, the S matrix is expressed in 
terms of the Hamiltonian matrix of the cavity, whose ele- 
ments are random from the unitary ensemble . From 
the resulting S and dS/ dE we compute the distribution of 
We observe that the probability density of $ is a nar- 
row peak centered around zero. The variance decreases 
quickly with N (see the right inset). Strikingly, the fluc- 
tuations are in good agreement (within the numerical 
error) with the analytical prediction (|10(1 . These results 
demonstrate unambiguously that the fluctuations of the 
nonlinear conductance of a chaotic cavity are not sym- 
metric under field reversal and that they depend strongly 
on the fluctuations of the screening potential. 

In single-channel conductors, the asymmetry is large 
[var($) ~ 0.2(e 3 /hS) 2 for N = 2, not shown in Fig. 
and of the same order as a linear conductance fluctua- 
tion. But as the number of channels increases the asym- 
metry rapidly becomes much smaller than a linear con- 
ductance fluctuation or weak localization correction due 
to the smallness of the DOS fluctuations [20J. Therefore, 
experiments on few-channel conductors are most promis- 
ing for the detection of the asymmetry described here. 



Conclusion. — We have investigated departures from 
the Onsager relations in mesoscopic systems in the 
nonlinear regime. Because of the screening potential, 
the weakly nonlinear conductance is asymmetric under 
magnetic-field reversal. We have determined the condi- 
tions under which such departures are experimentally ob- 
served. Our approach can be applied to systems that ex- 
hibit similar phenomena such as metallic nanowires |2lj 
and molecular junctions \22\ . 

This problem was suggested to us by B. Spivak and 
D.M. Zumbuhl. We acknowledge P.W. Brouwer, S. Pil- 
gram and P. Samuelsson for helpful discussions. This 
work was supported by the Swiss National Science Foun- 
dation, the EU RTN under Contract No. HPRN-CT- 
2000-00144, and the Spanish MECD. 

Note added. — During completion of this work we be- 
came aware of related work by B. Spivak and A. 
Zyuzin |2.'"{l| treating a diffusive cavity. 
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